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Abstract 

It is argued that in the hmit of extreme disorder AC hopping is dominated 

by "percolation paths". Modelhng a percolation path as a one-dimensional 

path with a sharp jump rate cut-off leads to an expression for the universal 

AC conductivity, that fits computer simulations in two and three dimensions 

better than the effective medium approximation. 
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While ordered solids show no frequency- dependence of their conductivity at frequen- 
cies below phonon frequencies, disordered solids are characterized by AC conductivity that 
varies as an approximate power-law of frequency The exponent is usually less than 

one, but often quite close to one. As the frequency goes to zero the conductivity stabilizes 
and becomes frequency-independent. The characteristic frequency marking the onset of DC 
conduction has roughly the same temperature-dependence as the DC conductivity. These 
features are observed universally for electronically conducting disordered solids like amor- 
phous semiconductors polymers 11,^, doped crystalline semiconductors at helium 
temperatures ||10| (where the random positions of the dopant atoms becomes important), or 



high temperature superconductors above |jTT|, as well as for ionically conducting disor- 
dered solids like glasses or polymers p|,H,||J^. 



The standard model for AC conduction in disordered solids is the hopping model IT4-II8 
The simplest version is hopping of non-interacting charge carriers on a regular lattice with 
random symmetric nearest neighbor jump rates; this is the model to be studied below. Al- 
ternatively, a macroscopic approach may be adopted by considering Maxwell's equations 
for a solid with spatially randomly varying conductivity ]T^|-pT[]. For both models the limit 
of extreme disorder may be studied by letting the temperature go to zero when the jump 
rates/the macroscopic conductivities are thermally activated with randomly varying activa- 
tion energies. It has recently been shown by computer simulations [p0| -|2^ that in the low 
temperature limit, the AC conductivity becomes universal in both models, i.e., independent 
of the activation energy probability distribution, p{E). The effective medium approxima- 
tion (EMA) for both models predicts the same universal AC conductivity in the extreme 
disorder limit (in more than one dimension). If a denotes the conductivity relative to the 
DC conductivity and s is a suitably scaled dimensionless imaginary frequency ("Laplace 



frequency"), the EMA universality equation [B0-E2|, first derived by Bryksin for the model 



of electrons tunnelling between randomly localized positions p!3| , |23[| , is 

a Ina = s. (1) 
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Computer simulations of the macroscopic model pO| , pT| showed that this equation works 
very well in two and three dimensions. The situation is different for hopping. While the 
existence of universality was confirmed in two dimensions, it was found that the onset 
of AC conduction is considerably smoother than predicted by Eq. (|I]). No simulations of AC 
hopping in three dimensions have so far been reported studying the extreme disorder limit. 

For hopping, the EMA is thus qualitatively correct by predicting universality in the 
extreme disorder limit, but quantitatively inaccurate. This is perhaps not surprising, 



since universality at extreme disorder is due to the dominance of percolation effects ||22 
and the EMA is a mean- field theory for percolation (which is a critical phenomenon). The 
EMA replaces the disordered solid by an "effective" homogeneous solid with characteristics 
determined by a self-consistency condition. Such an ordered medium cannot a priori be 
expected to accurately represent conduction along the optimal "percolation" paths of an 
extremely disordered medium p4| . 



An alternative to the EMA is the "percolation path approximation" (PPA) proposed for 
the macroscopic model [PT|. According to the PPA, in the extreme disorder limit the AC 



conductivity is equal to that of a one-dimensional model with a sharp activation energy 
cut-off; this explains the existence of universality The idea is that percolation paths. 



which at extreme disorder dominate conduction in more than one dimension, have two char- 
acteristics: They are essentially one-dimensional and they only involve activation energies 
up to the "percolation energy", E^, defined from the bond percolation threshold, by 

= rp{E)dE. (2) 

J —oo 



The percolation energy is the activation energy of the DC conductivity |2^j2^. For the 



macroscopic model the PPA leads to a = 5/ ln(l + 5), which is close to Eq. (|1|) and gives a 
good fit to simulations in more than one dimension ||21|| . 

The hopping version of PPA is not analytically solvable. Below we derive an approxi- 
mation to the hopping PPA utilizing the one-dimensional EMA, which is known to work 
well Despite the above objections against the EMA, this procedure does make sense. 



because in one dimension the "effective" homogeneous medium still is one-dimensional, of 
course. We show by computer simulations that the one-dimensional EMA, henceforth re- 
ferred identified with the PPA, gives a better representation of low-temperature AC hopping 
in two and three dimensions than Eq. (|1|). Finally, we briefly discuss the consequences of 
these findings. 

To arrive at the PPA, hopping in one dimension with a sharp energy barrier cut-off, 
p{E) = ioT E > Ec while p{Ec) > 0, is addressed. In the "rationalized" unit system 
where the conductivity for a homogeneous system is equal to the jump rate IQ, the EMA 
equation for the AC conductivity a{s) in one dimension [|l^,^,^^-|32| is [where s = iu is the 



Laplace frequency, F is the jump rate, and the brackets denote an average over the jump 
rate probability distribution] 

/ ^-^ \ =0. (3) 

The quantity sG is s times the diagonal element of the Green's function for a random walk 
on a one-dimensional lattice with uniform jump rate a [the "effective medium"]; sG is given 
ll,!^,!! by 



sG = [1 + ^] . (4) 



s 

We are only concerned here with relatively low frequencies where sG << 1. To lowest order 
in sG Eq. (^) may be rewritten [^ 

- = (^^) . (5) 

a \T + sGa/ ^ ^ 

The right hand side may be expanded as a power series in sGa, leading to 

i =E(-^Gar(F-("+^)). (6) 

n=0 

Since F = Fq exp(— where (3 is the inverse temperature, the average (f-("+i)) is easily 
evaluated in the low temperature limit: If (3 = l3/p{Ec), one finds to leading order in 1//? 
(f-("+i)) = F(Ec)~("+^V[(^ + 1)/^]- When this is substituted into Eq. (D the following 
equation is obtained 



Letting s go to zero we find ct(0) = (3 r{Ec). Introducing the dimensionless Laplace fre- 
quency, 



4 (t(0) 



(8) 

wlienever sG « 1 Eq. (|) implies jSsG = \Js/a, where a = a/a{0). Substituting this and 
T{Ec) = a{0)/P into Eq. (|^ finally leads to the PPA expression for hopping, 

\/d In 1 + y/ld = Vs . (9) 

Due to the factor in Eq. (H), as the temperature is lowered towards zero the condition 
sG << 1 is obeyed in a wider and wider range of frequencies around the onset of AC 
conduction. 

We have carried out computer simulations of low-temperature AC hopping in one, two 
and three dimensions using the Fogelholm algorithm [^] to reduce the AC Miller- Abrahams 
electrical equivalent circuit of hopping [p^-|36[] according to a recently proposed scheme 



To speed up the calculations the lowest jump rates were set to zero; it was carefully checked 
that this does not affect the conductivity. Figure one shows the results of computer simula- 
tions of low-temperature AC hopping in three dimensions. Results are shown for averages 
of 100 simulations of the AC conductivity at real Laplace frequencies for four different ac- 
tivation energy probability distributions at the following "reduced" inverse temperatures: 
(a) (3 = 80; (b) (3 = 160; and (c) (3 = 320. The full curve is the PPA (Eq. (g)) while the 
dashed curve is the EMA (Eq. (|l|)). Empirical rescalings of the data were allowed in order 
to focus only on the shape of the conductivity curves. Figure one shows that universality 
is approached as the temperature goes to zero and that the PPA gives a good fit to the 
universal AC conductivity in three dimensions. 

The universality may be studied without use of empirical rescalings by plotting the slope 
of the Log-Log plot, d Loga/d Logs, as function of a. This is done in Fig. 2 for data from 



computer simulations in one, two and three dimensions. The computer simulations in one 
dimension were carried out for systems with a sharp activation energy cut-off, to check the 
validity of Eq. (^). The full curve is the PPA prediction and the dashed curve is the EMA 
prediction in more than one dimension (Eq. (|1])). Clearly, the PPA works better than the 
EMA in two and three dimensions; in fact, the PPA works very well in three dimensions. 

As mentioned in the introduction, for the macroscopic model the EMA works very well 
in two and three dimensions [it is exact in one dimension]. For the macroscopic model the 
EMA universality prediction (Eq. (||)) is very close to the macroscopic PPA [^]. In view 



of the above presented results for hopping, it now appears that the EMA works well for the 
macroscopic model because the EMA happens to be close to the PPA, in contrast to what 
is the case for hopping. The two models have the common characteristic that, in the extreme 
disorder limit, conduction becomes essentially one-dimensional and that consequently the 
PPA gives a good description of the universal low-temperature AC conductivity for both 
models. 

To summarize, we have presented evidence showing that the PPA, despite being based on 
a naive one-dimensional picture of percolation, is a good model for AC hopping conduction 
in the extreme disorder limit in two and three dimensions. These findings have important 
experimental consequences. The low-temperature universal conductivity of the macroscopic 
model is different from that of the hopping model. The macroscopic model incorporates 
Coulomb interactions via Maxwell's equations, while hopping traditionally is concerned with 
non-interacting particles. Consequently, it is in principle possible to determine the relevance 
of Coulomb interactions by measuring the low-temperature AC conductivity. 
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FIGURES 



FIG. 1. 

Log-Log plot (base 10) of computer simulations (symbols) of low-temperature AC con- 
ductivity in three dimensions at real Laplace frequencies for four different activation energy 
probability distributions, compared to the PPA (Eq. (^), full curve) and the EMA (Eq. 

(H) , dashed curve). The dimensionless "reduced" Laplace frequency, s, is defined by Eq. 

(I) - however, an empirical rescaling was allowed to focus exclusively on the shape of the 
conductivity curve; a is defined by cr = a{s)/a{0). The jump rates are F = Fq exp(— /?£^), 
where the activation energy E is chosen randomly according to the following probability 
distributions ||21[: Asymmetric Gaussian [p{E) oc exp(— _E'^/2), < < oo] (x); Cauchy 
[p{E) oc 1/(1 + ^2), < E <oo] (+); Exponential [p{E) oc exp(-E), < E < oo] (O); and 
Box [p{E) = 1, < E < 1] (□). To speed up the calculations all jump rates with activation 
energy larger than E^ + 6.4//3 were set to zero [E^ is defined from the bond percolation 
threshold in Eq. (Q)]. In terms of the dimensionless inverse temperature f3 = f3/p{Ec) the 
figure shows data for 100 averages of simulations of cubic lattices with sidelength where 
(a) l3 = 80[N = 29]; (b) (3 = 160 [A^ = 54]; (c) (3 = 320 [A^ = 100]. 

FIG. 2. 
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The slope of the Log-Log plot of a{s) at real Laplace frequencies, d Logiocr/d Logios, 
as function of LogioO" for simulations (symbols, as in Fig. 1) at the inverse dimensionless 
temperature (3 = 320 in: (a) one dimension [100 averages of 8192 lattices]; (b) two dimen- 
sions [10 averages of 880X880 lattices]; (c) three dimensions [100 averages of 100X100X100 
lattices] for the four activation energy probability distributions of Fig. 1. The simulations 
in one dimension were carried out with a sharp activation energy cut-off at = 1 in order 
to show the validity of Eq. (^) in one dimension. In two and three dimensions, to speed up 
the calculations, all jump rates with activation energy larger than Ec + 6.4//? were set to 
zero. The simulations are compared to the predictions of the PPA (full curve, Eq. (^)) and 
the EMA (dashed curve, Eq. (|)). Both the EMA and the PPA predicts that the slope of 
the Log-Log plot goes to one as s — > oo, but the PPA works better than the EMA in two 
dimensions and much better than the EMA in three dimensions. 
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